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A  MATHEMATICAL  MODEL  FOR  PREDICTING  THE 
DAMPING  TIME  OF  A  MERCURY  DAMPER 

by 
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ABSTRACT.  The  mercury  damper  is  a  device  for  damping 
the  nutation  of  a  rotating  body.  This  report  explains  how 
a  mercury  damper  works,  and  derives  a  simple  expression 
for  the  damping  time. 
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FOREWORD 


Research  on  the  mechanism  of  a  mercury  damper  was  prompted  by  the 
need  for  a  nutation  damping  device  for  the  Air  Force  Lunar  Probe.  This 
report  is  the  result  of  a  theoretical  study  which  was  initiated  in  July 
1958  by  the  Systems  Analysis  Branch  of  the  Naval  Weapons  Center.  The 
study  was  carried  out  under  Advanced  Research  Projects  Agency  (ARPA) 
Project  Order  3-58,  and  completed  8  September  1958. 

In  this  analysis,  only  a  single  idealized  process  is  assumed  to 
occur.  Although  the  physical  situation  is  certainly  not  as  simple  as 
this,  measured  values  of  the  damping  time  have  compared  favorably  with 
those  predicted  by  the  theory.  More  experimental  work  is  needed,  however, 
before  the  predicted  mechanism  can  be  finally  accepted. 

This  report  is  released  at  the  working  level  for  informational 
purposes  only. 


Released  by  Under  authority  of 
C.  P.  SMITH,  Head  _  F.  H.  KNEMEYER,  Head 
Infrared  Systems  Division  Weapons  Development  Department 
20  November  1970 
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Wetted  area 

Resistance  coefficient  as  a  function  of  the  Reynolds  number 

Frictional  drag  force 

Component  of  H  along  P 

Component  of  H  along  T 

Total  angular  momentum 

Polar  moment  of  inertia 

Transverse  moment  of  inertia 

Length  of  channel  filled  with  mercury 

T  x  P 

Center  of  gravity 
Symmetry  axis 
Unit  vector  along  OP 
Reynolds  number 
Posi.ion  vector 
Component  of  R  along  N 

Component  of  R  along  P 

Component  of  R  along  T 

Unit  vector  in  plane  defined  by  OP  and  11,  ar.d  directed  so 

that  T  •  H  >  0  and  T  •  P  =  0 

Acceleration  in  direction  of  F  due  to  motion  of  coordinate 
system 


iv 


m  Mass  of  mercury 

r  Lever  arm 

t  Damping  time 

v  Peak  velocity 

Y  /ngle  between  H  and  P 

Yj  Final  value  of  y 

Yj  Initial  value  of  Y 

v  Kinematic  viscosity 
p  Specific  gravity 

to  Total  angular  velocity 
o)^  Component  of  to  along  P 
)p^  Initial  value  of  to^ 

u>t  Component  of  to  along  T 
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INTRODUCTION 


A  mercury  damper  is  a  device  consisting  of  a  circular  channel 
partially  filled  with  mercery  fixed  to  a  rotating  rigid  body  so  that 
if  any  nutation  of  the  rigid  body  occurs,  it  is  automatically  removed 
by  the  damper.  The  mercury  damper  accomplishes  this  by  extracting 
kinetic  energy  from  the  system.  The  time  required  to  stop  the  nutation 
is  the  damping  time,  which  is  characteristic  of  the  particular  combina¬ 
tion  of  damper  and  rigid  body. 

The  following  is  a  simplified  analysis  of  the  motion  of  the  rigid 
body,  and  of  the  effect  of  the  damper  upon  that  motion.  The  resulting 
equation  for  the  damping  time  will  aid  in  the  design  of  a  mercury 
damper  which  will  suit  the  requirements  of  the  designer. 


DESCRIPTION  OF  MOTION 


Consider  a  rigid  body  with  axial  symmetry  with  the  greatest  moment 
of  inertia  about  the  symmetry  axis,  spinning,  without  the  Influence  of 
external  forces  (see  Fig.  1).  Suppose  the  body  moves  in  the  most  gen¬ 
eral  way  so  that  the  symmetry  axis,  OP,  the  total  angular  momentum,  H, 
and  the  total  angular  velocity,  w,  are  non~coplanar. 

The  locatlon_of  a  point  in  the  rigid  body  shall  be  specified  by  a 
position  vector,  R,  drawn  from  the  center  of  gravity,  0,  as  an  inertial 
reference.  An  orthogonal  triad  is  defined  by  the  unit  vectors:  P  along 
OP,  T  in_the  plane  defined  by  0P_and_H*  ®nd  directed  so  that  T  •  H  >  0 
and  T  •  P  ■  0,  and  N  defined  by  T  x  P. 


R*RP  +  RN  +  R.T 
p  n  t 

and 

_  _  • 

w  •  w  P  +  Y  N  +  w  J* 

P  t 
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FIG.  1.  Spacial  Relationships  of  the 
Defined  Quantities. 
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EFFECTS  OF  DAMPER  UPON  RIGID  BODY  MOTION 


In  the  problem  that  is  being  considered,  a  coherent  mass  of  mercury 
is  allowed  to  flow  through  a  circular  channel  that  is  coaxial  with  the 
rigid  body  (see  Fig.  2).  The  frictional  drag,  F,  between  the  walls  of 
the  channel  and  the  mercury  produces  a  torque  about  OP  which,  by  Newton's 
second  law,  must  equal  Hp: 


ui  I 

P  P 


Letting  r  represent  the  lever  arm,  and  substituting  in  Eq.  1, 


Y 


rF 

H  sin  y 


(2) 


Note  that  if  the  damper  were  not  present,  then  Y  would  equal  zero,  and 
the  assumption  of  the  first  section  that  OP,  H,  and  u  are  non-coplanar 
would  not  be  true. 


FIG.  2.  Damper  as  Fixed  to  Rigid  Body. 


(Note  the  position  of  the  mercury  in  the  channel  relative 
to  the  total  angular  momentum  vector,  H.) 
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For  the  mercury  flow  rates  to  be  considered,  the  Reynolds  number, 
given  by  R  =  v1a>“1  (where  v  is  the  peak  velocity,  L  is  the  length  of 
channel  filled  with  mercury,  and  V  is  the  kinematic  viscosity)  is  suf¬ 
ficiently  large  so  that  the  flow  shall  be  considered  turbulent.  Under 
these  conditions,  the  frictional  drag  can  be  represented  approximately 
by  the  empirical  relation: 

F  =  ~  CfApv2  (3) 

where 

A  =  Wetted  area 
p  =  Specific  gravity 
v  =  Peak  velocity 

=  Resistance  coefficient  which  is  a  function  of  the 
Reynolds  number. 

For  this  problem,  Cf  lies  in  the  range  0.004  >  Cf  >  0.002.  The  v  is 
a  complicated  function  of  the  dimensions  of  the  channel  and  of  the 
position  within  the  channel.  For  this  reason,  v  is  assumed  to  be  equal 
to  the  average  speed  of  the  mercury  in  the  channel,  which  is  approxi¬ 
mately  true.  Now  that  an  expression  for  the  frictional  drag  is  available, 
the  damping  time  can  be  calculated. 


DERIVATION  OF  THE  DAMPING  TIME 


It  is  supposed  in  the  analysis  that  the  nutation  is  violent  enough 
so  that  the  mercury  is  forced  centrif ugally  to  a  position  in  the  channel 
which  is  far  from  H  (see  Fig.  2),  as  H  turns  about  OP  (or  OP  about  H) 
at  the  rate: 


co.H  I 

JL-R  =  _R 


co 


H.  Rouse  and  J.  Howe,  "Basic  Mechanics  of  Fluids,"  p.  170, 
New  York,  Wiley,  1953. 
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However,  the  rigid  body  turns  about  OP  at  Wp.  Therefore,  the  speed  of 
the  mercury  relative  to  the  channel  is 


-1  ,-l. 


v  =  r  (X  /It  -  1)ul  =  r(It  -I  )  H  cos  y 


Combining  Eq.  2  and  3,  and  substituting  for  v 


Y  =  - 


1  3  _  .  /T-1  -1.2  ,,2  2 

2  r  CfAp(lt  -  Ip  )  H  cos  y 

H  sin  y 


rearranging, 


Y  sin  Y  1  3„  A„,,/-r“l  ,-ly2 

2  "  "  2  r  CfApH(It  '  Xp  > 

cos  Y 


Integrating  from  Yi  to  Yf.  and  replacing  H  by  u)p^  Ip/cos  Yi  yields: 


1  - 


cos  Yj 
cos  Yx 


I  C£A0T3  Bpi  IpU;1 


c>2 


(4) 


This  is  the  time  required  for  the  frictional  drag  of  the  mercury  damper 
to  cause  the  rigid  body  to  precess  so  that  y  is  changed  from  Yi  to  yf. 
The  most  useful  case  is  the  one  for  which  Yf  =  0,  since  t  is  then  the 
damping  time.  However,  the  relation  is  not  strictly  true  when  Yf  is 
arbitrarily  small,  since  it  was  assumed  that  the  nutation  was  violent 
enough  to  constrain  the  mercury  to  a  fixed  position  in  the  coordinate 
system;  i.e.,  to  a  point  bearing  a  fixed  relation  to  H  and  OP.2  Under 
these  circumstances,  the  condition  that  F  =  ma  must  he  satisfied,  where 

1  2 
F  -  y  CfApv 


This  is  experimentally  observed. 
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is  the  frictional  drag,  m  is  the  mass  of  the  mercury,  and  a  is  the 
acceleration  in  the  direction  of  F  due  to  the  motion  of  the  coordinate 
system.  Wien  a  becomes  everywhere  smaller  than  F/m,  then  Eq.  4  is  no 
longer  valid. 

The  vector  acceleration  at  any  point  R  is: 


R=RP+RT  +  RN 
p  t  n 

From  Fig.  1,  neglecting  y, 


P  -  UtN  »  <Ip/It>  V 

f --<yy  “^--vv2  “pf 

N  =  (I/I.)  to  T  =  -(I  /IJ2  u>2  N 

p  t  p  p  t  p 

Substitution  in  R  yields 


R  =  Rn(I/IJ  <*>  a)  T  -  R.  (I  /IJ2  u>2  T  -  R  (I/IJ2  u)2  N 
pptpt  tpt  p  npt  p 


a  is  the  component  of  R  along  r  x  R,  where  r  *  RnN  +  RtT: 
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a 


(r  x  R)  •  (R) 
r  x  R 


-  (R  N  -  R  T)  •  (R) 
r  t  n  '  '  ' 


Multiplying, 


a  =  -  -  R  R  (I  /I)  ui  u) 
r  n  p  p  t  p  t 


I  to  =  H  cos  Y 
P  P 

■  H  sin  y 

t  t 


Rn  R  r2 

a  - 3 — 2 -  cos  y  sin  y 


For  a  given  value  of  Y»  the  largest  value  of  a  is  that  for  which 
Rn  =  -r. 

Substituting  for  Rj,,  and  setting  a  =  F/m  will  determine  the 
minimum  value  of  y  for  which  the  observed  damping  mechanism  is 
operative : 


2  2 

F/m  =  Rp  H  /It  cos  y  sin  y 

but 

1  2 
F  =  i  Cf Apv 


Substituting  for  v^,  and  solving  for  Y: 


Y  .  =  arc  tan 

min 


CfApr2(l  -  It/Ip)2 
2m  R 


(5) 
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In  most  problems,  it  :Ls  found  that  Ymin  sniall  enough  so  that 
Yf  of  Eq.  4  can  be  set  equal  to  zero  without  significant  error: 

1  -  cos 

damping  time  =  ^ ^ - 

-T-  C  Apr  u>  I  (I..  -  I  ) 

2  f  pi  p  t  p 

provided  Y  .  is  small 
min 

Figures  3  and  4  are  graphs  of  R  and  Cf  for  use  with  Eq.  5  and  6. 


(6) 
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(LOG  R) 


FIG.  4.  The  Resistance  Coefficient,  Cf(R). 
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